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On the partial differential equations of mathematical phjsics. 

By 

E. T. Whitt axes in Cambridge. 


§ 1 . 

Introduction. 


The object of this paper is the solution of Laplace’s potential equation 


d^v_ Q 

ox 3 oy z oz 7 


and of the general differential equation of wave-motions 


dx s dy s ^ 02 1 at* 7 


and of other equations derived from these. 

In § 2 , the general solution of the potential equation is found. 

In § 3, a number of results are deduced from this, chiefly relating 
to particular solutions of the equation, and expansions of the general 
solution in terms of them. 

In § 4, the general solution of the differential equation of wave- 
motions is given. 

In § 5, a number of deductions from this general solution is given, 
including a theorem to the effect that any solution of this equation 
be compounded from simple uniform plane waves, and an undulatory 
explanation of the propagation of gravitation. 


§ 2. 

The general solution of the potential eqnation. 

We shall first consider the equation 

d 21 + d ' v -o 
dx 1 ^ V dz* “ 

which was originally given by Laplace*). 


*) Memoirs sur la theorie de Vanneau de Satume , 1787 . 
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This equation is satisfied by the potential of any distribution 0 f 
matter which attracts according to the Newtoman Law. We shafl first 
obtain a general form for potential-fonctions and then shall ^ew that 
this form constitutes the general solution of Laplaces equation. From 

the identity Jjf 

! i r __ *a _ 

y ft.-.)*+‘O'" 1 ' “ ‘ 

we see that the potential at any point (n, y, r) of a particle of maaa m, 
situated at the point (fl* c )> i 3 




(2u 


m /*_______— r .— 

2* / (z+ixcoBu+iy sin tf) — (c-f-i a cos u-|-i b ainu) 

o 

which, considered as a function of x, y, z, is an expression of the type 

tn 

J* f^z+ix cos u-\~iy sin u, ii)du, 

o 

where f denotes some function of the two arguments 

z -f ix cos u -r iy sin « and w - 

It follows that the potential of any number of particles m 1 , m*, w t , 
situated at the points (fli&iO, ( a 2 ^s c 2 )> * * *> a A c *)> * s 411 ex 

pression of the type 

2rt 

f [f^z+ix cos u+iy sin n, w) + cos uJ r l V sijl w ) 

V + fk( zJ r ix cos tt +*2f sin u ) 1 

or 

Srt 

J 9 f(z+ix cos it+iy sin u, u)dUj 
0 

where f is a new function of the two arguments 

z 4- ix cos tt 4* iy sin u and u . 

In this way we see that the potential of any distribution of matter 
which attracts according to the Newtonian Law can be represented by an 
expression of the type 

2 it 

J*f{z-\-ix cos u+iy sin u, u)du . 


B - 3 
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The question now naturally suggests itself^ whether the most general 
solution of Laplace’s equation can be represented by an expression of this 
type. We shall shew that the answer to this is in the affirmative. 

For let V(x, y, 2 ) be any solution (single-valued or many-valued) of 


the equation 


3*7, 3*7, 3*7 n 


Let (x 9 , y 0J r 0 ) be some point at which some branch of the function 
V(Xj y, 2 ) is regular. Then if we write 


Xq -4 Z, 


Vo + % 2 — z Q + Z 


it follows that for all points situated within a finite domain surrounding 
the point ( j ? 0 , y 0 , j 0 ), this br an c h of the function Y(x, y, 2) can be ex¬ 
panded in an absolutely and uniformly convergent series of the form 

V a. a 9 + Oj X. + Y 4* Cj Z 4- a 5 X. 2 + 3T* 4* % Z 2 -f- d* YZ . 

+ e i ZX + f z XT+< h X' + *-: 


Substituting this expansion in Laplace s equation, which can be written 


o'V , d*v , c*v 
dX' t ar 1 * oz 5 


0, 


and equating to zero the coefficients of the various powers of Z, X and Z 7 
we obtain an infinite number of linear relations, namely 


a» 4* & 2 + Cm — 0, etc. 
between the constants in the expansion. 

There are -!■»(»—1) of these relations between the 4-(n-J-l) (n-i-2) 
coefficients of the terms of any degree n in the expansion of 7; so 
that only {(» + 2) — -~n(n—1)} or (2»+l) of the coefficients 

of terms of degree n in the expansion of Y are really independent. It 
follows that the terms of degree n in V must be a linear combination 
of (2n-f-l) linearly independent particular solutions of Laplace’s equation, 
which are of degree n in X, Y } Z. 

To find these solutions, consider the expansion of the quantity 

(Z+iX cos u~\-i Y sin u)" 

as a sum of sines and cosines of multiples of w, in the form 

(Z+tX cos u+i Yam u) n - g 0 (X, Y f Z) + ^(Z, Y, Z) cos u 

4- g 9 (X, Y, Z) cos 2u H-f- g n (X f Y, Z) cos nu 

4- /ij(Z, Y t Z) sin u -f \{X, Y } Z) sin 2« 4- 

+ h n (X, Y y Z) sin nu. 
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Now g (X, 7, Z) and h m (X, 7,Z) are together characterised by the fact 
that the highest power of Z contained in them is Z n ~ m \ moreover 
gjx, Y, Z) is an even function of 7, whereas h m (X, 7, Z) is an odd 
function of 7\ and hence the (2»+l) quantities 

g a ( X, 7, Z), g t (X, 7, Z), ■ ■ ■, K(X, 7, Z) 

are linearly independent of each other; and they are clearly homogeneous 
polynomials of degree n in X, 7, Z-, and each of them satisfies Laplace’s 
equation, since the quantity [Z+iX cos M+iTsin«)" does so. They 
may therefore be taken as the (2»+l) linearly independent solutions of 

degree n of Laplace’s equation. 

Now since by Fourier’s Theorem we have the relations 

\X 

g n {X, 7, Z) — -i- J [Z+iX cos tt+iTsin »)" cos mudu, 

0 

2 ?* 

h m (X, Y, Z) — ^ f [Z+iX cos «+iTsin »)" sin mudu, 

0 

it follows that each of these (2»+l) solutions can be expressed in the form 

/« z+iz cos u+i Y sin w, u) du 

o 

and therefore any linear combination of these (2»+l) solutions can be 
expressed in this form. That is, the terms of any degree n in the ex¬ 
pansion of V can be expressed in this form; and therefore V itself can 
be expressed in the form 

2tx 

0 

or 

j*F(z-^r cos u^riy sin u — 2 0 —ix 0 cosu— iy 9 sinw, u)du f 

o 

or 

2 * 

J*f(z+ix cos u~riy sin u, u) du y 

o 

since the r 0 + ix d cos u + iy sin u can be absorbed into the second 
argument w. 

Now V was taken to be any solution of Laplace’s equation, with no 
restriction beyond the assumption that some branch of it was at some 
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point a regular function — an assumption which is always tacitly made 
in the solution, of differential equations 5 and thus we Have the result, that 
the general solution of Laplace’s equation 


is 


v . d* r , a a v 
+ W + ^ T ' mm 




ix cos u+iy s in u 1 u) du, 


where f is an arbitrary function of the two arguments 

z + ix cos u + iy sin u and u . 

Moreover, it is clear from the proof that no generality is lost by supposing 
that f is a periodic function of u. 

This Theorem is the three-dimensional analogue of the theorem that 
the general solution of the equation 


is 


d 21ju d lZ o 

ox 1 + dy* 


7 — f(a+iy) + g(x—iy). 


§ 3. 

Deductions from the Theorem of §2; Particular Solutions; 
Expansions of the General Solution. 


1 °. Interpretation of the solution. We may give to the general solution 
just obtained a concrete interpretation, as follows. 

Since a definite integral can be regarded as the limit of a sum, we 
can regard V as the sum of an infinite number of terms, each of the type 

V r — f r (z+ix cos u r +iy sin u r ) 


each term corresponding to some value of u r . 
But this term is a solution of the equation 


where 


3 s V d * V 

dx}^Tz* 

X r =- x cos u r + y sin i* r , 
Y r =— — x sin u r -f y cos u r , 
Z r — s, 


so that (X rJ Y r} Z r ) represent coordinates derived from (x, y, z) by a 
rotation of the axes through an angle u r round the axis of z. 
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Tina we see that the general solution of Laplace's equation can It 
regarded as the sum of an infinite number of elementary constituents 7„ 
each constituent being the solution of an equation 

d r V r PJr _ o 
3X? ^ d# ’ 

and tie axes (X,, T r> Z T ) being derived from tie axes (x, y,s) by a 

simple rotation round tie axis of z. 

2" The particular solutions in terms of Legendre functions. It is 

interesting to see iow tie well-known particular solutions of Laplace's 
equationin terms of Legendre Sanctions can be obtained as a case of 

the solution given in § 2. ^ 

The particular solutions in question are of the iorm 

(cos 9) cos mcp and r* P™ (cos 9) sin m<p 

(n 0 , 1 7 2 , • * *, oo j m * 0 7 1 7 2 , * • * ? n) ; 

where (r, 9, <p ) are tie polar coordinates corresponding to the rectangular 

coordinates (x, y, s), and where 

.. (— 1 ')"* 3 in m 9 <r + m (sin-" 8 ) 

P »( C0S 9 ) - 77T~ d(coseT+” ‘ 

Now tie function P“ (cos 8) can be expressed by tie integral 

2 Jt 

m 

p m ^ _ (n-rM) — l)--(n + l) / ^ cos £-j_£ sin 9 coscos rmh 

0 

and thus we have 

^p r ( a »g)c 0 s» 9 >- (n+ ” )( " + ”~ i) "' ( " —/* i*+iYFW)<™*y 

cos Tnijf cos Tftcp dt? 


2 * 

_ (w+i»)(n+«t-!)—(«+1) iy^( 2 + iy x f + y z cos #) n cos /»(^— <p) 

0 

2 * 

_ (n^rm)(n-\-m — l) \y j*(z + ix cos u 4* iy sin u) n cos mu dix • 

0 

We see therefore that the solution r*P™ (cos 9) cos mtp xs a numerical 
multiple of 

In 

J'{z+ ix cos u + iy sin w)" cos mudu. 


B -7 
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Similarly the solution r»P“ (cos 8) sin m<p is a numerical multiple of 


C(jg -j- ix cos u t" iy sin u) n sin w# du. 

o 

From this it is clear that in order to express any solution 

Is 

jf(s + ix cos u + iy sin u, u) du 
0 

of Laplace’s equation, as a series of harmonic terms of the form 
r»P™ (cos 9) cos m<p and r*P” (cos 9) sin mtp, 
a is only necessary to expand the function f as a Taylor series with respect 
to the first argument z + ix cos « + iy sin u, and as a Former senes with 

respect to the second argument u. . , 

As an example of this procedure, we shall suppose it required to 

find the potential of a prolate spheroid in the form 

2n 

/«* -j- ix cos u + iy sin u, u) du , 

0 

and to expand this potential as a series of harmonics. Let 


*’+r | z * _a 

" " + e* u 


be the equation of the surface of the spheroid; and suppose that it is a 
homogeneous attracting body of mass M. To find its potential, we can 
make use of the theorem that the potential at external points is the same 

as that of a rod joining the foci, of line-density 5 that “» 1 

4 (C —- <1 ) 

2 it Vc 1 —a* _ , 

3iif Cl C (c 2 — a 2 — £*)<*£ 

8 J J *Z ?+ ‘ XC0S “ + ‘ y8intt 


IS 




-UL, fie-*-s’, m + d “- 

8 »(c— a)**' l ' 

x 0 

where B is written for s + ix cos u + iy sin u. 

Expanding the integrand in ascending powers of % , we have the 

potential in the form 

In 


3 M. C f_!_ 

2 7tJ 11 • 3 ■ 


B 


+ 3 • 5 • B s ^ 5 • 7 • 5* ~ 


■| du . 
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2 it 

_i_ r iit 
J * +1 


P^(CO80) 

/•+ 1 


this gives the required expansion of the potential of the spheroid in Legendre 

functions , namely the senes 

1 , fc»-a»)JP.(co»g) (c 1 — a 1 )' &) q_ . . .1. 

3 -^{r _ 3?' 1 —rr? ^ 5-7-r 1 / 

This result may be extended to the case of the potential of an ellipsoid 
with three unequal axes, by using a formula for the potential of an 

ellipsoid given by Laguerre*) . . -p , 

3°. 2he particular solutions of Laplace’s equation, which involve Bessel 

functions. We shall next shew how the well-known particular .solutions 

of Laplace’s equation in terms of Bessel functions can be obtained as a 

case of the general solution. The particular solutions in question are 

of the form _ ... T n N . 

^ z J m (Jcq) cos mcp and er~J m (jCQ) sm mcp 7 

where h and m are constants, and 2 , g, (p are the cylindrical co-ordinates 
corresponding to the rectangular co-ordinates x, y, z, so that 

x — q cos (pj y * Q si 11 *P • 

Now if in the solution 

et s J m (]cQ) cos vup 
we replace by ya l n0 


n. 

J m (hq) — ifcos (md — lq sin d) dd, 


we find after a few simple transformations that 

2 * 


jti r n,~\ «. cos mudu. 

«* a J m (hq) cos »»9 — - 57 - / ^ 


The other solutions which involve sin m<p, can be similarly expressed: 
we see therefore that the solutions 

^ s J m (lcQ) cos mp and e^ s J m (kQ) am mcp 


*) C. R., 1878. 


B-9 
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art numerical multiples of 

2 n 
% 

gi(j-tM*coiu + iy«iiw«) coa mudli 


and 


I e 


m 

/•“ 


e*(*+.•*«»«+»>»»»*») gin mu du 
respectively . It follows from this that in order to express any solution 


In 

S« 


z + ix cos u.+ iy sin u, «) du 


of Laplace's equation as a sum of terms of the form 

e kx J m Qt q) cos m 9 and ef* J m (h g) sin m <p, 

it is only necessary to expand the function f in terms of exponentials of its 
first argument z -f- ix cos u 4 - iy sin u, and as a Fourier series ivith respect 
to its second argument u. 

As an example of the use which may be made of these results, we 
shall suppose it required to express the potential-function 

F— 1 + er a J 9 (ff) + e~ 2 *J 0 (2g) + er u J 9 (Sq) + * * - 

(where z is supposed positive) as a series- of harmonic terms of the type 
involving Legendre functions: and also to find a distribution of attracting 
matter of which this in the potential. This can be done in the following 
way. We have 

F— 1 4 - e- - 7 0 (p) + <r 5s 7 0 (2g) + <r 3, c7 0 (3$>) 4 - 


2 3t 


JL J *J \ _J_ g-*-»*eo*«g-2(j + txcos«+ »yua«) J du 


2 n 


1 f* _ du _ 

2« J 1 (f + izcoiu + iyaiau) * 


But if t be any variable different from zero, and such that 
we have 



-T + T-^F+^T.-*»«£ + •••> 


t\ < 2xj 


where B lf B 2 , are Bernoulli's numbers. Therefore, so long as z is positive 
and \z 4“ ix cos u 4 - iy sin u\ < 2ir i. e., so long as z is positive and 
x * 4" 2T 4- z* < 4?r 3 we have 
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in 


F*»— I I-i- : —:-b 4-+ -It (^-{-ixcosM +iy sintt) 4-) iu 

2 nj l js —J— xx coat* 4“ l V aiat4 2 2! ^ J J 


or 


z+ ix coa-u 4* l V 

1 - T> /„™ <V\ « >r» / /i\ . ^ 


+ 27 ^Pi(cos^)- r 3 P,(c03<9>+ T jtr i P 3 (coad) + *.. 

and $is is required expansion ofVasa series of harmonics involving 
Legendre functions . 

Next, since 


we Have 


_i_i+i + y_ i +_i 

. ff -« 2 r *4"2wi* * — 2 niit 7 

nm I 


In 


V-± fdu[± + —. - i-rr : - +y\ 

2jt / L2 5 -{- txcoa w -j“ *y am 14 l 

c/ **1 


+ 


z 4- ix coa i* 4- iy sin u 4- 2nijr 

_i_n 

z 4- ix coa u 4“ iy ain t* — 2nt« ] J 7 


or 


f-^4- —- 1 —- 4- y 1 

2 l/.-r* —U f/* -i— 2 " JammJ 


_ 1 i 1 1 

yx s 4-y 2 4“^ s r \V X * 4“ 2 T 4“ (3+2ni«y s ]/x z 4“ y 5 4" ( 2 — 2 ni«y*j 7 


l 


and therefore 7 can 5e regarded as the potential due to a set of attracting 
masses placed at equal imaginary intervals 2 in along the axis of z. 


§ 4 

The differential equation tX 4- \X 4- %X — hr XX ' 

ox- oy oz - o t 1 

We shall next consider the general differential equation of wave- 
motions. 

aw aw aw , 2 aw 
a^ 2 ay s ‘ dz i * at* ? 

where h is a constant. 

Writing kt for f, this becomes 

aw aw aw aw 

dx z ‘ dy 1 ‘ Bz z ““ dt* 7 

which we shall take for the present as the standard form of the equation. 

In order to find the general solution of this equation, we follow a 
procedure analogous to that of § 2. Let V{x , y 7 z 7 t) be any solution 
(single-valued or many-valued) of the equation; and let (x 0 , y Q , z 0J 1 0 ) be a 
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-lace at which some branch of the function V is regular. Then if we 

Jntos-* + X,y-* +5 *-* + *, *-* + *,» P °/ 3lble 

to expand this branch of the function 7 as a power-senes of the lorm 

V-a 9 +a 1 X+b 1 7+c l Z+d 1 T+a t X i + b t T*+c t Z t +d t T'-+e i X7 

+ f t X2 + g. XT + ft, YZ + *, YT + l, ZT + a, X s + • • •, 

which will be absolutely and uniformly convergent for a certain finite 
frnmm of values of X, 7, Z, T. Substituting this expansion in the 
differential equation, which may be written 

3‘ 7 3* V 3*7 3*7 

dx 1 + ar* + dz* “ Tt*’ 

and equating to zero the coefficients of various powers of X, 7 and Z, 
we obtain an infinite number of linear relations, namely 

4“ ^j ““ dtc.j 

between the constants in the expansion. There are -|-(»—1) «(»+l) of 
these relations between the -g- (w+1) (« + 2) ( w +3) coefficients of terms 
of any degree n in the expansion of 7; 30 that only 

j {(n + 1) (n+2) («+3) - (n-l)n(»+l)} 

or 

(n 4" l) 5 

of the coefficients of terms of degree n in the expansion of V are really 
independent. It follows that the terms of degree n in F must be a linear 
combination of (n-r l) 2 linearly independent particular solutions of degree 
n in X, Y ? Z 7 T. 

To find these solutions, consider the expansion of the quantity 
(X sin u cos v 4* Y sin u sin v 4“ Z cos u 4" X)* 1 . 

If we first take the expansion in the form 

% 4“ 9% cos v 4* 9* cos 4-4- 9n coa nv 

4- \ sin v 4- sin 2v 4- -\-h n sm nv 7 

we have seen in § 2 that g 0J g 1} • * *, g n) * * *, A n , are linearly indepen¬ 
dent functions of X, Y, Z and T. Moreover, g m and h m are of the form 
sin" 1 u>< a polynomial of degree ( n — wi) in cos u, and therefore each of 
them contains (n-wfl) independent polynomials in X, Y, Z, T. Thus 
the total number of independent polynomials in X, Y 7 Z 7 T } in the 
expansion of 

(X sin u cos v 4- Y sin u sin v 4- Z cos u 4- T) n 
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in sines and cosines of multiples of u and v , is 

(»+l) + 2n H- 2(«— 1) + 2(n—2) 4-[-2 

or 

( n +l) * 

Now each of these polynomials must satisfy the equation 

d*v , d*V , d*V d'V 
dx x + dy* dz x “ 1 

since the quantity 

(X sin tt cos i? 4* X" sin w sin v 4- Z cos u + T) n 

does so: and therefore they may be taken as the (n+1)* linearly inde- 
pendent solutions of the equation 

d * Y _l ill — d * V 

dx 3 dy 3 "* Si 1 

which are homogeneous of degree » is X, 7 ; Z, T. Now by Fourier s 
theorem we have 
2* 

* - ifa sin m cos v 4- 7 sin u sin v + Z cos u 4- T) n cos mv dv\ 


and since g m is of the form 


71—771 

^ u r sin ,tt u cos r u, 


where t* is one of the polynomials in question, it is clear that g m can 
be expressed as a sum of sines or cosines of multiples of u, according 
as m is even or odd; and the coefficient of one of these sines or cosines, 
say of cos sw, is 


if 9 - cos su du. 


It follows that each of the polynomials u r can be expressed in the form 


it 

f L/w du > 

0 

where f(u) denotes some periodic function of w; that is, it can be ex¬ 
pressed in the form 

2 it it 

sin u cos v 4- Y sin u sin v 4- Z cos u 4- T) n f(u ) cos mv du dv . 
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It follows from this that each of the (w-{-l) 2 polynomial solutions of 
degree n can be expressed in the form 

In n 

JI (X sin u cos v 4“ Y sin it sin v -}- Z cos u Yy* f(u, t?) du dv } 


o a 


where f(u , v) denotes some periodic function of u and t?; and therefore 
the terms of degree n in V can be expressed in this form. 

The function V itself can therefore be expressed in the form 


In n 


f ff(* sin u cos t? + [Tain u sin v + Z cos u + T } u f v ) dudv , 


«/ t/ 

O 0 


where f denotes some function of the three arguments 

X sin u cos v Y sin u sin v 4- Z cos u ■+■ T, u y and u; 

and f may without loss of generality be supposed to be periodic in u and v. 
Now 

X sin u cos u + Y sin u sin v + % cos u ■+■ Y 

— (x sin u cos v + y sin w sin i; H- cos u -J- f) 

— (a? 0 sin it cos v + y Q sin u sin v -j- cos u + Q] 

and the termo 

(x 0 sin it cos v -f y Q sin u sin v -J- r 0 cos u 4- f 0 ) 

can be absorbed into the arguments u and t?; moreover V was taken to 
be any solution of the partial differential equation; we have, therefore, 

on writing y for t 7 the result that the general solution of the partial 

differential equation of wave-motions 7 

d*V , d'V , 2*V ,, 

dx 1 + dy 1 + ds 2 “ K *** 9 


is 


it 

r -J Jf( x sin u cos v -j- y sin u sin v + z cos u + -p, w, vj du dv, 

o o 

where f is an arbitrary function of the three arguments 

x sin u cos v + V sin u sin v + z cos u + -t-, u an d v. 


alath«ra»tisch« Ar. n aian L V i_L 


23 


B- 14 
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§ 5. 

Deductions from the general solution of § 4. 

1°. The analysis of wave-motions. We shall now deduce from the 
general solution thus obtained a result relating to the analysis of those 
phenomena which are represented by solutions of the equation 

d*Y will. 

dx 1 dy 1 dz 1 dt 1 

If we revert to the fundamental idea of the definite integral as the limit 
of a sum of an infinite number of terms, we see that the general solution 


2 * it 


r-fff(sm u cos v + y sin u sin v + z cos u ■+■ -j-, u, vjdudv 


o o 

can be interpreted as me aning that V is the sum of an infinite number 
of terms of the type 


f(x sin u cos v + y sin u sin v + s cos u + ; 


u. 



there being one of these terms corresponding to every direction in space 
given by the direction-cosines 

sin u cos v, sin u sin v, cos u. 

The solution F can therefore be regarded as the sum of constituent 
solutions, each of the type 


F (x sin u cos v -f y sin u sin v -j- z 


cos u 



where the function F varies from one direction •(«, v) to another. 

Now let us fix our attention on one of these constituent solutions F. 
If for some range of values of the quantity ^ 

x sin u cos v y smu smv z cos u + -j*, 

the function F is finite and continuous, we can for this range of values 
express F by Fourier’s integral formula in the form 

as b 

dXJ F(ot) cos | X (x sin u cos v + y sin u sin v -j- z cos u -j- — Aa j dX da , 

0 a 



where a and 5 are the terminals of this range of values; or supposing 
the integration with respect to a to be performed, 


J aw™ {>■ (x sin u cos v + y sin u sin v -f- z cos u + ^ ? 


where g(X) denotes some function of X. 
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Now let os again revert to the idea of the definite integral as the 
limit of a sum. Then this latter integral can he regarded as the sum 
of aa infinite number of terms of the type 


™{x(x a mucosv + ysmu3mv + z<:osu + T )}, 


each term being multiplied by some factor depending on A. 

The solution V can therefore be regarded as constituted by the super¬ 
position of terms of this last type. But a term of this type represents 
a simple uniform plane wave; for on transforming the axes so that the 
new a-riw of x is the line whose direction-cosines are 


the term becomes 


sin u cos v, sin u sin v } cos u, 


cos 

sin 


(*+t)> 


which represents a simple plane wave whose direction of propagation is 
the new axis of x. We see therefore that (he general finite solution of the 
differential equation of wave-motions, 


ill ill 4. r ' v 

ox 1 dy 1 " r sF 


,, o-r 


can be analysed into simple plane waves f represented by terms of the type 

w, v ) ^ | A (x sin it cos v + y sin u sin v + z cos u + p*) } ■ 

It is interesting to observe that Dr. Johnstone Stoney in 1397*) 
shewed by physical reaso ning , and without any reference to the equation 

d*v , d*v , d'T _ /2 d s V 

dx 1 dy 1 dz 2 *“ at 2 

that all the disturbances of the luminiferous ether arising from sources 
of certain kinds can be resolved into trains of plane wares. 

2°. Solution of the equation 


d'V d*V , d*V v 
dx £ + dy * + dz 2 + K 


0. 


If a solution W of the equation 


a*T r fl w 

Jx* + dy 1 + dz 2 *“ dt 1 

be of the form FV* ; where F is a function of x, y, z only, which does 
not involve t , then V clearly satisfies the equation 

dx 2 + dy 2 + dz 2 ^ y ’ 


*) Philosoph . Magazine, (V) X LU L 
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and therefore, on reference to the general solution of the way e-motion 
equation found in §4, we see that the general solution of the equation 


w+w + S+ r -° 


18 


it 


T—JJ***>•*) f(u, v)dudv 
0 0 


3°. Deduction of the known particular solutions of the equation 

0. 


3*7 , 3*73*7^- 

i i i y 


dy 1 ' 3z* 

It is known that particular solutions of the equation 

3^7 3^7 3fJ 

3x* + 3r + 3* a + K U 


exist, which are of the form 

V-r-±J m+± (r) P;(co 3 0)“ s m<p 
2 

(w =»0,1,2,• • *; w —0,1,2,•••,»), 

where r, d, cp are the polar coordinates corresponding to x, y y z. We shall 
now shew how these may be derived from the general solution of the 
equation which has just been found. 

For let the general solution be written in the form 


lit X 


7 -J ’ f( u> gin u dudV) 


where f(u 7 v) is an arbitiary function of the two arguments u and v , 
which may without loss of generality be taken to be periodic in u and v. 

Now let the function f(u , t?) be expanded in surface-harmonics of u 
and v, so that 


2jt it 

, / / gi(- rain « co, * + y* i a 

n » 0 


u am 9 + a cos u 


> Y n (ti, v) sin u dtt dv 


where Y n is a surface-harmonic of order n, L e v if 

5 q sin u cos v, r\ =— q sin it sin v, £ =■ p cos w, 
are regarded as the co-ordinates of a point in apace, then g n Y n (u, v) is 
a homogeneous polynomial of degree n in §, rj, £, satisfying Laplace's equation 

a *7 ,3*7 3*7 

+ 3^ + 3? 7 * 0 * 
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Next, let the variables be changed by the substitution 

COB u — cos S cos o 4" sin Q sin <d cos v, 

sin u sin ([<p — v) » sin o sin v , 

sin » cos ( 9 - 0 ) - cos o sin 9 - sin a cos v cos 0, 

so that (a sin <d cos v, g sin o sin v, g cos o) are the co-ordinates of the 
point ( 5 , v, 5 ) referred to new axes, the line whose direction-cosines are 
(sin 9 cos 9 , sin 0 sin 9 , cos 6) being taken as the new axis of z. 

Thus 

2 it it 


-iff' 


gi'rcM® Y n (Uj v) sin a dm dv\ 


But a surface-harmonic of any order n remains a surface-harmonic of 
order » under any transformation of axes in which the origin is unchanged: 
and therefore TJu, v) is a surface liarmonic of order n in o and v ; 
and consequently it can be expanded in the form 

A.(9, v )P,(emz,) + A l J9, V ) PJ(cosci)cos®' + Al(d,<p)Pl(c o««)cw2®' 

-j-h A" (9, tp) P"(cos<o) cos nv 

+ S' (9, <p) P; (cos «) sin *' + ••• + K («, <P)K^° S <a ) shinv '’ 
where A m (6, 9 ), • • •, <p) are functions of 9 and 9 . Substituting this 

value for YJu, v) in "the integral, and performing the integration with 

respect to v y we have 

• ✓** 

V—^A,(9, 9 ) / e ir ' M "P„(cosaj) sin oiai; 

nmQ U 

0 

and in virtue of the relation*) 

i iV t (r) 

( -j\ 2 a +-*- 

•y) --; 

this can be written in the form 

F - Y r_T J »+4 ( r ) 9>) 

a = 0 

where /* m ( 0 7 <p) denotes some function of 0 and cp . 

*) A proof of this and several related results mil be found in a paper shortly 
to be published by the author. 
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Since the surface-harmonics Y n (d, cp ) were independent of each, other, 
the functions f n (9, cp ), will be independent of each olher and therefore 
each of the quantities 


r * <p) 

a 


will be a solution of the equation 

JP + JF + JF + V ~ 0 ’ 

But on transforming this equation to polar co-ordinates, and substituting 
the expression 

r ~ Tj n+ 1 (r)W, <?) 

for V, we find that the function f n (B, 9 ) must satisfy the differential 
equation for a surface-harmonic in d and 9 of order «. It follows that 
f n (6 , <p) can be expanded in the form 

f,(9, <P ) — A, p, (cos 9) + Al cos cpP* (cos 9) -)-V A n cos nip P“ (cos 9) 

+ B\ sin <p P* (cos 9) H-(- B" smn<pP”(cos 8), 

and thus the particular solutions 

r'T J 1 ( r) p;(cos&) mcp 


* + T 


are obtained. 

Moreover, it is clear from the above proof that in order to expand 
any solution 


J it it 


F- J J, ^(**to.co..+,.u.u.in.+ 3in u du dl] 

0 0 


of the equation 

as a series of the form 


f!I + i!I + £!I + 7^o 
a* 1 + a«* + oi‘ + v u 


'2 r ~ t 

71 = 0 2 


wfore IT* is a surface-harmonic of order n in 6 and 9 , i£ is owfy necessary 
to expand the function f(u 7 v) in surface-harmonics of u and v. 

4°. Expression of the solution of the equation 

g + % + Q + v.o 

as a series of generalised Bessel functions. 
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Another analysis of the solutions of the equation 

*r±£Z±zz.+v-o, 

IF + 3 y* + 3 *’ + 

entirely different from that given in 3°, can be found in the following way. 
Consider tlie expression 

,M-4)K)-t4-v)HH’K). 

if this expression be regarded as a function of a and t, it can for 
values of s and t be expanded as a senes of (positi 
negative) integral powers of s and t, the coefficients in this senes ben* 
functions of x^and r. Let the coefficient of the term m s”t" be deno 
by J^fx, y, s): so that we have the relation 

(*4) (‘4) -T* (*4) (' 4) + T*('4) y, *) • 

oa Tt™-*® 


This equation can be regarded as a generalisation of the equation 

.t-R) 

Jl» — <» 

which defines the ordinary Bessel functions; and we shall consequently 
can the functions J^(x, y,») generalised Bessel functions. 

We now proceed to establish some properties of the functions 
J (x, y, *); it will be seen that they are very similar to those of the 

ordinary Bessel functions. 

In the first place, since the expression 

r _.HR)RHdR)KK-H) 

satisfies the equation 

IF + dy‘ + 32* + 

it follows that each of the functions J^.O, V, *) soM«s equation 

iZ + m + ZI + v-o. 

ax' + dy 1 + 3* 

In the second place, we shall obtain an expression torJ^&y,*) as a definite 
integral By Laurent’s theorem, we know that the coefficient of s^ in 

the expansion of 

>H)HK’R)RK-K) 


T 3 .. "*0 
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ji, J' s - n -i e T z (‘-r) ( ,+ t) (*-t) + t'(' + t) dSj 


where C is any simple contour in the 5-plane surrounding the origin 5 
and again applying Laurent's theorem ; the coefficient of 2 " in this ex¬ 
pression is seen to be 


*//■ 
C D 




-i e 7*(-T) ( <+ t)-T»(*-t)( ,+ t) + 7’(' + t) 


where D is any simple contour in the 2 -plane surrounding the origin. 
Now write s =— e*“, t — e**. Thus we Aave 2 Ae result 

Sit Sit 

hSS- — min — + ixiinuco*f+ tyainutin*-!- i jct>» u 


0 0 


ichich may be regarded as the analogue of BessePs integral 

X 

™ “ J* 30S (nu — z sin u) du. 

0 

The functions J’ IJhn (x,y, 2 ) likewise possess an additiontheorem: for we have 

g l(* + «)(,-l) (1+4) -4^+*) (—7) ('"t) + T (,+<r) ( J+ T) 

s t'(-t) ( <+ t) ( # -t) 


* e 

x e 


and so 


,W-t)KK-B)HM-K) 

+ y + b, * + c)s m r 

m, ji 

“ ^ J m,n (*, 

Ji ntf it 

Equating coefficients on both sides of this equation, tee have the result 
+a,y+b ,2 + c) ^ (*> y> 2 ) Jm - r ,«)» 

=* — as £=* — o» 

w/iicA is 2 /ie addition-theorem for the generalised Bessel functions , and is 
the analogue of the well-known result 

£ « 
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We dull now shew tow the generalised Bessel functions furnish an 
analysis of the general solution of the equation 


3»7 

ax 1 




ay 


0 . 


For the general solution is, bj 2°, 


I* 


J v)dudv : 


where f(u, v) can without loss of generality he taken to he a periodic 

^^ow^et the function f(u, v) he expanded hy the extended form of 
Fourier’s theorem, in the form 




Then we have 


m*-® *» —' 


ft In 


^7 ( f gi(2smMCOi® + y»ia“»*“» + ,co, “ + m “ 

y “2j Zi *’*'*/ / 

0 o 


+**) dudv . 


Comparing this with the form just found for the generalised Bessel 
functions, we see that the general solution of the equation 

d'v . £!! + £IT + f-o 

3 l*' r + az> + 

can Jc 

A 

7-2 ^ 9 >*)’ 


jn»-« 


where the quantities «... are arbitrary constants. This furnishes an 

alternative analvsis of the solution to that given in 2 . 

5®. Gravitation and Electrostatic Attraction explained as modes of Wave- 

disturbance. . , 

The result of 1°, namely that any solution of the equation 

3*7 d ' y j- I 

can he analysed into simple plane waves, throws a new lightion the 
nature of those forces, such as gravitation and electrostatic attrac ion, 
which vary as the inverse square of the distance. For if a system of 
forces of this character he considered, their potential (or their componen 
in any given direction) satisfies the equation 


T 3 - 
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B*V , 

gar 3 ‘ gy 3 g * 1 


0 . 


and therefore d fortiori it satisfies the equation 

d*V d*V £7 ,,^7 

g ® 3 ' gy 3 g ^ 3 gtf 1 

where Tc is any constant It follows from 1° that this potential (or force- 
component) can he analysed into simple plane waves in various directions, 
each wave being propagated with constant velocity. These waves interfere 
with each other in such a way that, when the action has once been set 
up, the disturbance at any point does not vary with the time, and 
depends only on the coordinates ( x , y, z) of the point. 

It is not difficult to construct, synthetically, systems of coexistent 
simple waves, having this property that the total disturbance at any point 
(due to the sum of all the waves) varies from point to point, but does 
not vary with the time. A simple example of such a system in the 
following. 

Suppose that a particle is emitting spherical waves, such that the 
disturbance at a distance r from the origin, at time t f due to those waves 

whose wave-length lies between — and , ■, is represented by 

[L (t t uj)> 


2 dii sin (pVt — fir) 
xp r 


where V is the velocity of propagation of the waves. Then after the 
waves have reached the point r, so that (Yt — r) is positive, the total 
disturbance at this point (due to the sum of all the waves) is 



Take pVt — pr 


or, since 


it is 


— y, where y is a new variable. Then this disturbance is 



r 


The total disturbance at any point, due to this system of waves , is therefore 
independent 'of the time , and is everywhere proportional to the gravitational 
potential due to the particle at the point . 
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It is clear from t3ie foregoing tliat the field of force due to a 
gravitating body can be analysed, by a “spectrum analysis” as it were, 
into an infinite number of constituent fields; and although the whole field 
of force does not vary with the time, yet each of the constituent fields is 
of an undulatory character, consisting of a simple wave-disturbance propa¬ 
gated with uniform velocity. This analysis of the field into constituent 

fields can most easily be accomplished by analysing the potential — of 
each attracting particle into terms of the type 

sin (ft Vt — pr) 
r 

as in the example already given. To each of these terms will correspond 
one of the constituent fields. In each of these constituent fields the 
potential will be constant along each wave-front, and consequently the 
gravitational force in each constituent field will be perpendicular to the 
wave-front, L e. the ivctves iviU be longitudinal 

But these results assimilate the propagation of gravity to that of 
light: for the undulatory phenomena just described, in which the varying 
vector is a gravitational force perpendicular to the wave-front, may be 
compared with the undulatory phenomena made familiar by the electro¬ 
magnetic theory of light, in which the varying vectors consist of electric 
and magnetic forces parallel to the wave-front. The waves are in other 
respects exactly similar, and it seems probable that an identical property 
of the medium ensures their transmission through space. 

This undulatory theory of gravity would require that gravity should 
be propagated with a finite velocity, which however need not be the same 
as that of light, and may be enormously greater. 

Of course, this investigation does not explain the cause of gravity; 
all that is done is to shew that in order to account for the propagation 
across space of forces which vary as the inverse square of the distance, 
we have only to suppose that the medium is capable of tr an s m itting, with 
a definite though large velocity, simple periodic undulatory disturbances, 
similar to those whose propagation by the medium constitutes, according 
to the electromagnetic theory, the transmission of light. 




